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ABSTRACT 



Molecular line polarization is a unique source of information about the magnetic fields and anisotropics in the circumstellar envelopes 
of evolved stars. Here we present the first detection of thermal CO(J = 2— >1) and SiO(7 = 5^4, v = 0) polarization, in the envelope 
of the asymptotic giant branch star IK Tau. The observed polarization direction does not match predictions for circumstellar envelope 
polarization induced only by an anisotropic radiation field. Assuming that the polarization is purely due to the Goldreich-Kylafis 
effect, the linear polarization direction is defined by the magnetic field as even the small Zeeman splitting of CO and SiO dominates 
the molecular collisional and spontaneous emission rates. The polarization was mapped using the Submillimeter Array (SMA) and is 
predominantly north-south. There is close agreement between the CO and SiO observations, even though the CO polarization arises 
in the circumstellar envelope at ~ 800 AU and the SiO polarization at < 250 AU. If the polarization indeed traces the magnetic field, 
we can thus conclude that it maintains a large-scale structure throughout the circumstellar envelope. We propose that the magnetic 
field, oriented either east-west or north-south is responsible for the east-west elongation of the CO distribution and asymmetries in 
the dust envelope. In the future, the Atacama Large Millimeter/submillimeter Array will be able to map the magnetic field using CO 
polarization for a large number of evolved stars. 
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1. Introduction 

Molecular line polarization observations can provide invaluable 
insights into the magnetic field and/or anisotropics of the cir- 
cumstellar environment of asymptotic giant branch (AGB) stars. 
^ If even a relatively weak magnetic field is present, the molecular 
. gas in the envelope of the AGB star will show linear polarization 
■ when the magnetic sublevels of the rotational states are exposed 
to anisotropic emission (the Goldreich-K ylafis effec t , Gold reich 
& Kylafis 1981, 1982). Alternatively. iMorris et all (Il985h pre- 
dict polarized line emission arising from molecules in an enve- 
, lope that have a preferred rotation axis because of radial infrared 
emission from the central star. In this case, deviations from ra- 
dial symmetry in the polarization angles could be due to a non- 
spherically symmetric envelope and/or to a (non-radial) mag- 
; netic field. While polarization observations of molecular lines, 
in particu lar of CO, iii star-forming regions are r elatively com- 
mon (e.g. ICortes et al.ll2005l; iBeuther et al.ll2010 '). the only pre- 
viously published significant detection of molecular line polar- 
ization around an A GB star is that of th e CS(7 = 2— > 1) transition 
around IRC+10216 (iGlenn et al.lll997h . 

Most information about the magnetic fields in the envelopes 
of AGB stars is currently derived from ma ser observations (e.g. 
IVlemmings et al.lf2005l : iHerpin et al.ll2006h . The maser observa- 
tions reve al strong magnetic field s throughout the entire enve- 
lope (e.g. Vlemmings et al. 7011). As maser observations typi- 
cally probe only a limited number of lines-of-sight through the 
stellar envelope, the polarization of thermal molecular lines can 
more easily reveal the circumstellar magnetic field morphology. 



A strong magnetic field is, in addition to both the binary and disk 
interactions, one of the possible mechanisms for shaping the typ- 
ically spherical AGB winds into an aspheiical planetary nebula 
(e.g. Bahck & Frank 2002; Frank et al. 2007|). 

Here we present Submillimeter Array (SMA) polarization 
observations of the molecular Unes in the circumstellar enve- 
lope of IK Tau, which is a well-studied M- type AGB star. It has 
a period of 500 days (Kukarkin et al.' 1971), and a fairly high 
mass-loss rate of 1x10"^ Moyr"' (Ramstedt et al. 2008). The 
star has a rich circumstellar chemistry, as indicated by the de- 
tection of a large number ( >20) of different molecular species 
(see e.g., iKim et aljQoiOt iDecin et alJl2010L for the latest re- 
sults). The circumstellar gas distribution has been mapped in 
CO(y = 1^0), CO(y ^ 2--»l) (Castro-Can-izo et al.l 12010). and 
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CO(y = 3^2) (iKim et alJ [20101) emission. ICastro-Carrizo et al.l 
(1201 Oh found a western elongation in its innermost region (of 
~ 1000 AU), and on larger scales a flattened circumstellar enve- 
lope that is elongated in the east-west plane over ~ 10.000 AU. 
Elongated, elliptical circumstellar structures are also found on 
scales from a few to a few te ns of AU in maps of the SiO (e.g., 
ICotton et al]l2010l) and H2O (iBains et alJl2003h maser emission. 
The SiO maser observations indicate that there i s a magnetic 
field close to the central star of 1.9-6.0 Gauss dHerpin et al.l 
.2006.) . 

2. Observations and data reduction 

The observations of IK Tau were done with the SMA on 
October 8 2010 in the compact configuration with the lower- 
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and upper-side band (LSB and USB) covering the frequency 
ranges of 216.8-220.7 GHz and 228.9-232.8 GHz respectively. 
This frequency setting was specifically selected to cover the 
CO(7 = 2— >1) transition at 230.538 GHz, but also covered tran- 
sitions of SiO, SiS, and SO. We list the detected molecular tran- 
sitions in Table [1] The correlator setup provides a spectral res- 
olution of ~ 0.8 MHz, which corresponds to ~ 1.0 km s ' at 
230 GHz. At 230 GHz, the beam size is ~ 3.0 x 2.6 arcseconds 
with a position angle of 8 1 .4°. The observations lasted for a total 
track of nine hours and included observations of the phase cali- 
brator J0423-013, bandpass and polarization calibrator 3C454.3, 
and the primary flux calibrators Callisto and Neptune. Based on 
the flux measured for J0423-013 (~ 2.2 Jy beam"'), we esti- 
mate that the total flux uncertainty is < 10%. The right ascen- 
sion (a) and declination (S) of the phase center for the obser- 
vations of IK Tau were taken to be o-jiooo - 03''53"'28'.84 and 
5j20oo = -(-11°24'22".56. The data were calibrated initially in 
the MIRIAD softwar e package using the sp ecific polarization 
calibration technique (iMarrone & Raoll2008l) . Using a standard 
SMA polarization observing schedule, the polarization calibra- 
tor 3C454.3 was observed regularly, covering > 120° in paral- 
lactic angle. The gain calibrator J0423-013 is linearly polarized, 
but does not show any circular polarization. The derived ampli- 
tude gain corrections thus do not affect the amplitude calibration 
of the right- and left-circular polarizations produced by the SMA 
observation. The antenna leakage solutions were found to be up 
to 2% in the USB, near the CO(7 - 2^1) line and up to 5% in the 
LSB, near the SiO(y - 5— >4, v = 0) line. The linear polarization. 
Pi, and electric vector polarization angle (EVPA) determined on 
3C454.3 (LSB: P, = 0.71 + 0.12%, EVPA= -47 ± 6°; USB: 
Pi = 0.80 + 0.09, EVPA= -52 + 2°) indicate that after calibra- 
tion the remaining leakage is < 0.5%. After calibration of the 
flux, phases, bandpass and polarization, we exported the data 
into the Common Astronomy Software Applications (CASA) 
package and performed self-calibration on the strongest spec- 
tral channel of the SiO line in the LSB. Although the SiO emis- 
sion is only marginally resolved, we used the initial total in- 
tensity image as a model. This assumes that the SiO emission 
has no significant circular polarization, which was verified to 
be the case. The SiO self-calibration solutions were compared 
with self-calibration solutions derived from the strongest spec- 
tral channel of the CO line (also following an initial image as a 
starting model) in the USB and were found to be identical. We 
then applied the SiO self-calibration solutions to both LSB and 
USB as the SiO solutions were obtained using data with higher 
signal-to-noise ratios. Maps of the total intensity Stokes I and 
the linear polarization intensities Stokes Q and U were created 
using 'briggs' weighting. To optimize the polarization detection, 
spectral channels were averaged over an interval of 5 km s"' for 
the CO, SiS, and SO lines and 2.5 km s ' for the SiO fines. 
As a results, the channel root-mean-square (rms) noise ranges 
from ~ 35 mJy beam ' for the CO line to ~ 65 mJy beam ' 
for the SiO line. We note that the noise in the LSB is ~ 30% 
higher than that in the USB. We also imaged the CO line with- 
out self-calibration in order to test the effect of self-calibration 
on the polarization. The EVPAs and linear polarization structure 
were found to be the same as when self-calibration was applied. 
Only the significance of the polarization detections decreased by 
~ 15%. 

3. Results 

Emission was detected from the CO(7 = 2— >1), SiO(7 = 5— >4, 
V = 0), SO(y = 6^5), and SiS(y = 12^1 1) transitions, hereafter 



Table 1. Properties of observed lines 



Line 


vo 


Fr 


F, 


pmax 




[GHz] 


[Jy beam ' ] 


[Jy beam"' 
kms-'] 


[%] 


SiO(y = 5^4, v = 0) 


217.105 


15.2 


379 


8.8 


SiS(y= 12^11) 


217.818 


3.7 


93 


<5.4 


SO(J = 6^5) 


219.949 


2.0 


58 


<6.2 


CO(J = 2-»l) 


230.538 


7.5 


204 


13.1 



indicated by CO(2-l), SiO(5-4), SO(6-5), and SiS(12-ll). The 
spectra of these transitions is shown in Fig. [T](left). Significant 
polarization was detected for the CO(2-l) and SiO(5-4) lines 
shown in Fig. |2] and online Fig. [3] To determine the significance 
of the weak linear polarization signal, we carefully removed 
the positive bias, which was introduced when the Stokes Q and 
U measurements were combined to produce the linearly polar- 
ized e mission (Pi - -\JQ^ + U^), following I Wardle & Kronberg 
('1974'). As the errors in Pi are not Gaussian-distributed, the 
significance is defined by determining the probability intervals. 
The contours in Fig. |2]and online Fig. |3]thus correspond to the 
3cr - 99.73% and higher probability intervals. The peak of the 
linear polarization corresponds to 4.9cr and 6.8cr detections for 
the CO(2-l) and SiO(5-4) lines, respectively. The maximum po- 
larization fractions are listed in Table{T] The EVPAs for both 
CO(2-l) and SiO(5-4) are predominantly oriented north-south. 
The polarization angle of CO(2-l) ranges from approximately 
-10° in the southern part of the envelope in the channel around 
Vlsr = 25 km s ' to ~ 10° in the northern part in the chan- 
nels with Vlsr - 25 and 35 km s"'. A somewhat smaller spread 
(-7 to 4°) is seen for the SiO(5-4) line, although for this line the 
polarized emission is confined close to the peak of the Stokes I 
emission. The measurement uncertainties for the EVPA ranges 
from 6-10° for CO(2-l) and 4-10° for SiO(5-4). This indicates 
that the difference in EVPA between the northern and southern 
part of the CO envelope is potentially real. 



4. Discussion 

4.1. Polarized molecular line emission 

The models of polarized molecular line-emission prescribe that 
the polarization can be caused by a preferred radiation direc- 
tion (Morris et al. 1985) or an interaction with a magnetic field 
(Goldreich & Kylafis 1981, 1982). In the first case, the stel- 
lar radiation will lead to polarization predominantly perpendic- 
ular to the radial direction. When the magnetic field is non- 
radial, and strong enough to determine the molecular alignment 
axis, competition between the field an d radiation direction will 
change the polarization characteristics. Morris et al.](ll985h state 
that for any field stronger than a few juG, the magnetic field in 
the envelope determines the molecular alignment axis. For the 
Goldreich-Kylafis effect, the linear polarization is oriented along 
the magnetic field as long as the Zeeman splitting dominates the 
collisional and spontaneous emission rates. 

The importance of the magnetic field in determining the po- 
larization clwacteristics of molecular line emission was esti- 
mated in iKvlafisI (11983.) . We used their Eqs. 1, 2, and 3, with 
the parameters adjusted to describe the CO(2-l) transition, and 
assumed a magnetic field strength ~ 1 fiG in the CO(2-l) region, 
a magnetic moment fi = -0.27 jun (with fi^ the nuclear magne- 
ton), and a reduced dipole matrix element d = 0.11 debye. It 
is then clear that a hydrogen number density < 3 x 10^ cm"-^ 
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Fig. 1. (left) Spectra of the observed molecular lines in the envelope of IK Tau. The emission is summed over all SMA baselines. The 
vertical line indicates the Vlsr - 34 km s~' stellar velocity, (right) The stellar continuum and dust emission (color) and integrated 
CO(2-l) line emission (contours) of IK Tau. The contours are drawn at 10, 20, 40, 80, and 160 Jy beam ' km s ' . 





Fig. 2. The polarization of the CO(2-l) line at 230.538 GHz in the circumstellar envelope of IK Tau. Channels are averaged across 
intervals of 5 km s ' width and are labeled according to the velocity at the lower end of the spectral bin. The color indicates the CO 
emission and the contours are the hnearly polarized intensity. Contours are drawn at a significance of 3 and 4cr (debiased, see text). 
The line segments indicate the electric vector polarization angle (EVPA). The beam size is indicated in the lower left corner of each 
panel. 



is sufficient for the magnetic field to determine the direction of 
the polarized CO(2-l) emission. As the SiO(5-4) emission orig- 
inates from closer to the star, the stronger magnetic field will 
also result in a dominant Zeeman splitting. Owing to the field 
strength measured around I K Tau and other AGB star s (> 1 mG 
within several hundred AU, IVlemmings et al.l (l2005h X the lin- 
ear polarization thus likely traces the magnetic field morphology 
when assuming the polarization originates fro m the Goldreich- 
Kylafis eff^ect. As indicated in Kylafij (Il983h the CO EVPAs, 
arising from the Goldreich-Kylafis eff'ect, are determined by the 
magnetic field direction on the plane of the sky. This is shown 
irrespective of the direction of the velocity gradient in the enve- 
lope, except that the relative optical depth caused by the velocity 
gradient determines whether the EVPAs are parallel or perpen- 
dicular to the magnetic field direction. 



The observed linear polarization fraction is, however, 
larger than p redicted from the general Goldreich-Kylafis effect. 
ICortes et all (^005) predict < 3% for CO(2-l) when the opti- 
cal depth is low (t ~ 1). They also show that higher fractional 
polarization can be achieved when a strong anisotropic radia- 
tion source is present. However, in that case, the polarization 
direction is determined by the anisotropic radiation field that, 
in the case of a circumstellar envelope, would cause tangen- 
tial polarization. As the linear polarization in the envelope of 
IK Tau is neither tangent ial nor purely radial as predicted by 
iDeguchi & WatsonI (1 19841) from multi-level calculations involv- 
ing infra-red excitation, it is unlikely that the anisotropic stel- 
lar radiation field contributes significantly. The high polarization 
fraction is potentially caused by other anistotropies in the cir- 
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cumstellar enviroment. More detailed modeling will be required 
to determine the origin of the large fractional polarization. 

4.2. The magnetic field of IK Tau 

For both the CO(2-l) and SiO(5-4) lines, we observe polariza- 
tion that is slightly offset from the peak of the total intensity 
emission. For CO(2-l), the largest polarization fraction is found 
in the south part of the blue-shifted side of the envelope and in 
the north part of the red-shifted side. The observed structure in 
the polarized emission is likely an effect of similar structure in 
the optical depth, as the largest fractional polarization is found 
when the optical depth t ~ 1 . 

As the circumstellar magnetic field strength has been shown 
to be sufficient to determine the molecular alignment axis, the 
polarization vectors are either parallel or perpendicula r to the 
magnetic field direction in the framework outlined by KvlafisI 
( 119831) . In that case, the overall field geometry is predomi- 
nantly either east-west or north-south. As only relatively few 
independent polarization vectors are measured, a more accurate 
magnetic field reconstruction will require more sensitive obser- 
vations. We are able to compare the EVPA direction derived 
from the CO(2-l) and SiO(5-4) EVPAs with the pr evious mea- 
surem ents made of the SiO maser polarization by Herpin et al.l 
(l2006l) . In that paper, the 86 GHz SiO maser linear polarization, 
averaged over several individual maser features with the 30-m 
IRAM telescope, was found to have an EVPA of 170-180° for 
the masers blue-shifted with respect to the star The EVPA then 
rotates to 140-160° for the red-shifted masers. Thus, the CO(2-l) 
and SiO(5-4) EVPAs are, for the blue-shifted emission, remark- 
ably consistent with the avera ged 86 GHz S iO masers. However, 



for the 43 GHz SiO masers, iCotton et al. (2010) found, from 
high-angular resolution interferometric observations a much 
more complex and variable polarization structure. The CO(2- 
1) polarized emission is detected about 3" (~ 800 AU) from 
the star and that of SiO(5-4) originates within ~ 1"(~ 250 AU), 
which c ould be an indication tha t, as in the case of the supergiant 
VX Sgr d Vlemmings et aLluOl lb . the magnetic field has a large- 
scale component that is preserved throughout the envelope. As 
the field is then mainly oriented either east-west or north-south, 
it might be directly related to th e slight east-west extent see n in 
the circumstellar CO envelope (ICastro-Carrizo et all 1201 0[ and 
Fig. [pright)) and the ea st- west asymmetry in the dust distribu- 
tion dWeiner et al.ll2006 l). Magneto-hydrodynamical simulations 
indicate, for example, that a dipole magnetic field in a circum- 
stella r envelope can re sult in equatorially enhanced density pro- 
files dMatt et alj|2000l) . However, we note that the elliptical SiO 
mas er structure is orien ted more in a northeast-southwest direc- 
tion dCottonetalJl2010l) . 



4.3. Prospects for ALMA 

Observations such as the ones presented here will undergo a 
revolution when the Atacama Large Millimeter/submillimeter 
Array (ALMA) becomes fully operational. With 50 antennas and 
in typical weather conditions, ALMA will reach the same sen- 
sitivity in the linear polarization as reached with the SMA (as- 
suming similar angular resolution and spectral binning) in less 
than two seconds. In particular for IK Tau, it will be able to 
improve the observations to a lOcr linear polarization detection 
with higher spatial resolution (1" x 1") and velocity resolution 
(2 km s"') in less than 20 min on-source observing time. While 
the SMA observations filter out a large fraction of the CO(2-l) 



flux, ALMA with the compact array and total power antennas 
will recover most of the emission and reach a similar polariza- 
tion detection level even for the evolved stars with weaker CO 
emission. Thus, on the basis of the SMA observations presented 
here, with ALMA it will be possible to map out the CO(2-l) and 
other molecular-line linear polarization in great detail for a large 
number of (post-)AGB stars and planetary nebulae in a relatively 
short time. As a consequence, ALMA polarization observations 
will provide unique information on the magnetic field morphol- 
ogy around evolved stars. 

5. Conclusions 

We have presented the first observations of polarized 
CO(7 = 2^1) and SiO(7 = 5^4, v = 0) emission in the circum- 
stellar envelope of an evolved star. As the polarization is nei- 
ther radial nor tangential, the linear polarization does not match 
typical predictions for linear polarization induced solely by an 
anisotropic radiation field or optical depth effects. It thus po- 
tentially traces the magnetic field even when influenced by 
anisotropic emission from the central AGB star IK Tau. In that 
case, the magnetic field is oriented either east-west or north- 
south and a good correspondence is found between the direction 
derived from the CO polarization at ~ 800 AU from the cen- 
tral star and that from the thermal SiO at ~ 250 AU. The slight 
east-west elongation of the CO and the previously observed 
dust asymmetry could be related to the large-scale magnetic 
field morphology. Further observations and theoretical work are 
needed to properly understand the influence of the magnetic field 
and the envelope and radiation anisotropics on the circumstellar 
linear polarization. It may then be possible to fully describe the 
circumstellar magnetic field structure. In particular, ALMA will 
be able to uniquely map out polarization for a significant sample 
of AGB and post- AGB objects. 
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Fig.3. As Fig|2]for the SiO(5-4) v = line at 217.105 GHz. In this case, channels are averaged over intervals of 2.5 km s ' and 
contours are drawn at debiased polarized intensity levels of 3, 4, 5, and 6cr. 



